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ABSTRACT 

Observational consequences of tidal disruption of stars (TDEs) by supermas- 
sive black holes (SMBHs) can enable us to discover quiescent SMBHs, constrain 
their mass function, study formation and evolution of transient accretion disks 
and jet formation. A couple of jetted TDEs have been recently claimed in hard 
X-rays, challenging jet models, previously applied to 7 -ray bursts and active 
galactic nuclei. It is therefore of paramount importance to increase the current 
sample. In this paper, we find that the best strategy is not to use up-coming X- 
ray instruments alone, which will yield between several (e-Rosita) and a couple of 
hundreds (Einstein Probe) events per year below redshift one. We rather claim 
that a more efficient TDE hunter will be the Square Kilometer Array (SKA) 
operating in survey mode at 1.4 GHz. It may detect up to several hundreds of 
events per year below 2 ~ 2.5 with a peak rate of a few tens per year at z ^ 0.5. 
Therefore, even if the jet production efficiency is not 100% as assumed here, the 
predicted rates should be large enough to allow for statistical studies. The char¬ 
acteristic TDE decay of however, is not seen in radio, whose flux is quite 

featureless. Identification therefore requires localization and prompt repointing 
by higher energy instruments. If radio candidates would be repointed within a 
day by future X-ray observatories (e.g. Athena and LOFT-like missions), it will 
be possible to detect up to ~ 400 X-ray counterparts, almost up to redshift 2. 
The shortcome is that only for redshift below 0.4 the trigger times will be less 
than 10 days from the explosion. In this regard the X-ray surveys are better 
suited to probe the beginning of the flare, and are therefore complementary to 
SKA. 
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Introduction 


Since the late 70s it has been snggested that stars torn apart by the grav itational he ld 
of a snpermassive b 


Frank fc Rees 


1976; 


ack hole (SMBH) may be observed as flares from Earth flHills 


Rees 


1988 


1975 


Phinnev 


1989I) . These are called tidal disrnption events 


(TDEs). These flares would be caused by sudden accretion of the star debris, which would 
feed the SMBH at an ever decreasing rate, M oc This theoretical expectation is for a 


com plete disruption of a star in parabolic orbi t , after at least several days 


(e.g. 


Lodato et ah 


2009 


Havasaki et ah 


2013 


Guillochon fc Ra.mirez-R,niz 


is expected to be independent on the ratio of pericenter to tidal radius fiSari et al 


Stone et ah 


2013|). 


rom the peak 


2013), and 1 


2010 


The detection and study of these flares can deliver important astrophysical information. 
On the one hand, they allow us to detect otherwise quiescent SMBHs and estimate their 
masses. This would inform theory of galaxy-SMBH cosmological co-evolution. On the 
other, they constitute a unique opportunity to study the - highly theoretically uncertain 
~ formation of an accretion disc and its continuous transition through different accretion 
states. As the accretion rate decreases, we can in principle observe a disc which transits 
from an initial super-Eddington phase, lasting several months, passing through a slim and 
later a thin disc regime, and ending its life, years later, in a radiative inefficient state. 
The super-Eddington phase -which occurs only for SMBH masses 10 ^ Mrr^ — is highly 


uncer tain, but it may be associated with a copious radiative driven wind flRossi fc Begelman 

-1 r- 


2009 

h which therma 

y emits ~ 10 ^^ — 10 


fStrnbbe & Onataert 

2009; 

Lodato & Rossi 

2011 ) 


2nil|). The disc lum inosity 10^^ — 10^° erg 


s ^), instead, peaks in far-UV/soft X-rays (ILodato fc Rossill2nill) . Of paramount theoretical 


importance would also be the possibility to investigate the formation and evolution of an 
associated jet, powered by this sudden accretion. There is no specihc theory for the jet 
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emission from TDEs. Astronomers mainlv assume a phenomenological description (e.g. 


Van Velzen et al 


2011 


7 -ray bursts (e.g. 


Canizzo et al 


Metzger et al 


2012 


20111 ) or borrow theory de veloped for blazars and/or 


Tchekhovskov et al. 


20131) . In general, non-thermal 


emission in X-rays and radio is the jet signature. 


Handful of Candida 
ROSAT all sky survey 


Survey flGezari et al 


(IVan Velzen et al. 


e TDEs 1 0 ) have been c 


flKomossa 


2009 


2002 : 


Donley et al 


Gezari et al. 


2012 


electe d so far, particularly in 


20021). in GAL 


Gampana et al, 


;X Deep Imaging 


2011 ) and in SDSS 


2011al) . These “soft” events are believed to be associated with the 


disc and wind thermal emission. The presence of a bright optical flare in the initial 
super-Eddington months makes optical surveys a useful tool for discovery. Signihcant 
advances in optical transient surveys are expected to be achieved by the Panoramic 
Survey Telescope and Response System (Pan-STARRS) and the Large Syno ptic Survey 


Telescope (LSST). Two can didates have been cla imed in Pan-STAR RS data flGezari et al. 


2012 


Ghornock et al 


(iHoloien et al. 


2014 ). three in PTF data f Arcavi et al 


20141) and one in ASAS-SN 


20141 ). but the total number expected seems to be much higher. For example 
in the Stt Survey, claims in literature range from 200 to ~ 1557, while in the medium deep 


survey there is more con sensus that ~ 15 — 20 should be found (jStrubbe fc Quataert 


Van Velzen et al 


2009 


2011a|) . Thousands of candidates could be, instead, detected by LSST, 


with its 6-band (0.3 — 1.1 micron) wide-held deep astronomical survey of over 20000 square 
degrees of the southern sky, usin g an 8.4-meter ground-based telescope flStrubbe fc Quataert 


2009 


Van Velzen et al. 


2011af) . However, these estimates are probably upper limits. 


because galactic nuclei can heavily absorb optical light. 


More recently, two candidates TDEs were triggered in the hard X-rav band bv ' 

le BAT 

instrument on board of Swift ( 

Bloom et al. 

2011 ; 

Burrows et al. 

2011 ; 

Genko et al. 

2012 

)• 


A multi-frequency follow-up from radio to y-rays revealed a new class of TDEs, where we 
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are likely observing the non-thermal emission from a relativistic jet. The jet emission is 
responsible for the hard X-ray sp ectrnm (with power-law slope (3 ~ 1.7) and the increasing 


radio activity flLevan et ah 


20111) . detected a few days after the trigger. 


Given the lack of statistics and of a solid theoretical framework for the non-thermal 
emission, we will take the best stndied of these two events, Swift J1644-I-57 (Sw J1644 
in short), as a prototype for the study presented in this paper, where we investigate the 
detection capability of both SKA and future X-ray observatories. 


Sw J1644 w as hosted by a s 


with its center (jZauderer et ah 


ar fo rming galaxy at z = 0.354 and in positional coincidence 


20111) . Its X-ray peak luminosity ~ 3 x 10^® erg s“^ was 
reached after a couple of days from the trigger, and it persisted at the level of > 10^^ erg 
s“^ for about 1 year. During its decay, the X-ray emission was approximately described by 
a temporal law, the same as that expected for the fallback of stellar debris (see Figure 
[T]). After ~ 500 day from the trigger, the X-ray flux declined by two orders of magnit ude 


and it has been associated with a shut off of the relativistic jet (IZanderer et ah 


20131). 


The modelling of the X-ray luminosity suggests that Sw J1644 is associated wit h a light 
supermassive black hole ^ IO^Mq (e.g. 


Burrows et al. 

2011; 

Ganizzo et al. 

2011) 


Variability at optical wavelengths within the host was not detected, while transient 
emission was seen in infrared, becoming stronger at longer wavelengths, especially at 
millimeter and radio wavelengths. Radio (1.4. and 4.8 GHz) observations from Westerbork 
Synthesis Radio Telescope (WSRT) showed a bright source. EVLA observations of the radio 


Lorentz factor T ~ 2 of the outflow ( 

Zauderer et a 

2011 

). The radio lightcurve disp 

ays a 

rebrightening starting one month after the trigger (J 

Sereer et al. 

2012; 

Zauderer et al. 

2013) 


The emission peaks around s everal months. 


at 600 days after the trigger fIZauderer et al 


bllowe d by a decline. Radio observations stop 


20131) . The radio behavior is not compatible 
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with the blast wave model borrowed from 7 -ray bursts by iMetzger et ahl (120121) . and 
indicate a m ore complicated jet s t ructu re, like perhaps in the magnetically arrested model 


proposed by 


Tchekhovskov et ah 


(120131) . Snapshot rates of jetted TDEs in radio band have 


been computed for the hrst time by Van Velzen et ah (2011). Differently from their work, 
we adopt here a different modelling for the radio lightcurve and a more detailed one for 
the black hole mass function, which includes the redshift dependence. We also account for 
a stellar mass function. We broaden up our investigation to include X-ray detection and 
follow-ups. 

Finally, a 200-s quasi-periodic oscillation (QPO) was detected by both Suzaku and 


XMM, ~ 10 and 19 days after the Swift/BAT trigger, respectively ( Reis et al.ll2012[) . QPOs 


are regularly detected in stellar mass BHs, but there is no hrm physical interpretation of 
these phenomena. However, most models strongly link the origin of high-frequency QPOs 
with orbits or resonances in the inner accretion disk close to the BH. This may cause 
variable energy injection into the jet, which consequently results in variability in the X-ray 
emission. This in terpretation led to estimate a BH mass between 5 x IO^Mq and 5 x 10 ®Mq 


(IReis et al 


201211 . 


In this paper, we predict the detection rate of jetted TDEs considering current and 
future radio surveys (NVSS -I- FIRST, VLT Stripe 82, ASKAP, VLASS and SKA) and 
X-ray instruments (Swift, eRosita, Einstein Probe, Athena, LOFT). In addition, we discuss 
the ability of these instruments to constrain important physical parameters. 


The paper is organized as follows. In §2, we take Swift J1644 as a prototype and we 
describe our phenomenological model for X-ray and Radio emissions. In §3, we discuss the 
black hole distribution functions used in this paper. In §4, we present our Monte Carlo 
calculations. Our rates for current and future surveys are presented in §5. A summary and 
implications of our results can be found in § 6 . Finally, we draw our conclusions in §7. 
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Throughout this paper we use the following cosmological parameters: = 0.25, 

= 0.75 and Hq = 70 km/s/Mpc. 


2. Modelling the Lightcurve 

A tidal disruption event of a star by a SMBH causes a transient accretion disc to 
form, whose accretion rate is set by the rate at which the stellar debris falls back to the 
black hole under its gravitational pull. How matter circularizes to form a disc and whether 
this process is accompanied by outflows and their characteristics are snbject to intense 
investigations, as mentioned above. From phenomenology and theory, we know that in the 
presence of an accretion disc and some ordered magnetic field, matter and energy outflows 
in form of (relativistic) jets are produced. In the absence of fully consistent simulations of 
jet production by a tidal disruption event, we use below a simplified description for the jet 
energy content as a function of time. This is partially supported by analytical and numerical 
calculations (see references above) and partly by the observed features of the X-ray emission 
of Sw J1644. In particular, its temporal decay (~ suggests that at least in this 

optically thin regime, the X-ray luminosity scales as the accretion rate. As a consequence, 
it supports a scenario in which the star was completely tidall y disrupted, since partial 


disru ption would lead to a shallower decay of the fallback rate flGuillochon fc Ramirez-R.uiz 


20131) . Moreover, a partial disruption is difficult to reconcile with a long lasting super 


Eddington accretion phase, which may be needed to power the jet for its total duration of 
~ 500 days. Finally, the modelling of the X-ray luminosity suggests that Sw J1644 is the 


hole 10 ’^Mq (e.g. 


Burrows et ah 

2011; 

Canizzo et ah 

2011) 











2.1. Jet kinetic power 


We work in the framework of two identical jets, with < 1/r. The total energy 
injected in the two jets is Lj = ejMfbC^, where Cj is the jet production efficiency, which we 
assume constant in time, and the gas fall back to form a disc occurs at a rate Mfb. For a 
complete disruption of a star in parabolic orbit the fallback rate can be approximated by 


Mfb('r) = iifp 


(Rees 

1988; 

Phinnev 

1989) 


^min + “T 


-5/3 


( 1 ) 


1989|). The lag time “r” is the time from the beginning of the debris 


accretion, that roughly happens after a time 


^min 


41 day. 


from the star disruption, in the galaxy rest frame. More precisely, fmin is the minimum 
time it takes the most bound debris to come back to pericenter after the star has 
been torn apart. Here and in the following, Mq is the BH mass in units of IO^Mq 
and the mass of the disrupted star in units of IMq. The peak of the accretion 
ratqj is quite intuitively the mass of the star divided by the characteristic timescale, 

Mp ~ (1/3) ~ 1-9 X g s“^. In our description, the jet is launched at 

the onset of accretion (f = 0), as there are no strong theoretical reasons why it should be 
delayed. The temporal evolution of the jet energy is thus 


where 


= h 


^min “I” ^ 


-5/3 


.2 _ 1 ^7 .. in45 _ ^-1 I ^ ^-l/2^1/2_ 


Lj^p = ejMpC ^ 1.7 X 10 erg s 


0.01 


( 2 ) 


(3) 


^In the formula used in this pap er, we assume the stan dard linear relation between mass 


and radius of the star. See eq.6 in (ILodato fc Rossi 


201ll). 
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Note that the larger the black hole mass, the lower the peak luminosity, because the 
characteristic timescale increases. Viceversa, the jet luminosity decreases with m*. 


2.2. X-ray 


The unabsorbed 1-10 keV lightcurve of Sw J1644 is shown in Fig|T] (black dries). 


Activity was already detected b y BAT 3 days be: 


beginning of XRT observations flBurrows et al 


2011 


ore the BAT “official” trigger and the 
Therefore there is an indication that 


the trigger (i.e. when the hrst photon was detected by XRT) happened approximately 
r Ri 3 day after the actual disc and jet formation. The observed time interval r is related to 
the rest frame analogous quantity hy t = f{l + z) and in this case r ~ 2 day. Accounting 
for this delay, the general behaviour of the X-ray lightcurve as a function of time At since 
the trigger (At = 0) can be reproduced by 


Lx,iso (At) PS 1.5 X 10"^® erg s 


-1 


r -|- At 


r 


-5/3 


(4) 


(Fig. [T], solid line). Specihcally, Lx,iso is an isotropic equivalent luminosity, computed 
from the X-ray flux. Note that here r = 3 is a hxed time delay, unlike r in eqsil] and [2l 
Superimposed to this baseline trend, there is a complex structure of flares and dips where 
the flux oscillates within two orders of magnitude in the hrst ten days of observations. It is 
clear that eg |4] does not capture this large variability , possibly associated with jet precession 


and nutation fjSaxton et al 


20121 : 


Stone et al. 


2OI2I) . But in absence of a compelling theory 


for these sudden X-ray variations, we prefer to reproduce the upper part of the envelope 
that contains the initial variability, since the BAT instrument was triggered by one of the 
peaks in the lightcurve. We will discuss later how this choice affects our X-ray TDE rate 
estimates. 


The Swift/XRT (0.3-10 keV) spectrum of Sw J1644-I-57 is well described by an absorbed 
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power-law with a photon index 13 ~ 1.7 — 1.8 and iVn ~ 2 x 10^^ cm“^ flBurrows et al.ll201ll) . 
The observed BAT spectrnm at early times and its count rate later on (up to the beginning 
of June) are consiste nt with an extrapolation at higher energies of the XRT spectrum 


(iBurrows et ah 


201lh . This suggests that we are observing the same component in both soft 


and hard X-ray bands. The average spectrum is consistently hard (1.4 ^ (3 1.7) during 


the whole emission, al 


initial variable phase fjSaxton et al. 


hough a spectral s oftening is observed during the short dips in the 


2 OI 2 I) . The radiation efficiency in 1-10 keV banc 


the fraction of the total luminosity emitted in that band) is ~ 0.20 flBurrows et al 


jTe 


201lh . 


With this last information, we can calculate the associated jet kinetic luminosity from the 
observed light curve, once we assume a jet opening angle 6-^ and a Doppler factor 5, 

L'i = Tx,iso(l - cosej)/{eJ^). 


With the highest probability, our line of sight is at an angle ~ T"^ (i.e. the inverse of 
the Lorentz factor T) that grazes the relativistic beam, and h T. The fact that there are 
no sharp breaks in the lightcurve may indicate that the whole emitting area was visible, 
i.e. T”^ > 6y Therefore, we further assume a jet opening angle of a similar size of the 
relativistic beaming, say 6j ~ r“^/2, and we get a jet power at the trigger time (At = 0) 
of Tj(r) ~ 1.5 X 10"^^ erg s“^(r/5)“'^. Since Tj(r) = ejMfb(r)c^, it turns out that to have 
an efficiency cj greater than 1% requires m* < IMq, for T < 5. In particular, m* = IMq 
gives efficiency between roughly 1% and 37% for 2 < T <5, that are in agreement with 


numer ical simulations of jets from highly super-Eddington accretion discs 


Sadowski et al. 


(120141) . Lower mass stars would give a higher efficiency range. We therefore assume in the 


following that Sw J1644 is the result of the disruption of a solar mass star. However, it is 
clear that this is just a tentative, though reasonable, choice, since the stellar mass cannot 
in fact be univocally determined, unless we can actually measure 6j. 


Assuming Sw J1644 as a prototype, we can adopt a general description of the X-ray 
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lightcurve in the 1-10 keV band, when we catch the flare after a time r from the beginning 
of the event, 


-hx,iso(^^) -hxjt 


r + At\-=-'= 


r 


(5) 


The (isotropic eqnivalent) Inminosity at the time of the trigger (At = 0) is 

— cos6*j) ~ Lj(r)ex2 (r/6*j)^, which can be written more explicitly as 


Lx,t 1.63 X 10"^® erg s ^Mg 


X 


^min + T 


-5/3 


( 6 ) 


where t = r 1(1 + z) and the radiation efficiency e:^{z) varies becanse of the spectral shifting 
with redshift, 

. (.) = 0 20 {E.(l + - (Edl + 

where we assnme (3 = 1.8, Ei = 1 keV, E 2 = 10 keV and Zgw = 0.35. We note that this 
correction is in the sonrce rest-frame and applies to nnabsorbed flnxes. 


In eqini we set 23.9. Indeed, any combination of these qnantities that gives a 

factor 24, allows us to reproduce the Sw J1644 X-ray luminosity at the trigger time. The 
degeneration should then be lifted, when we need to choose a Lorentz factor to compute the 
TDE rates. From the X-ray luminosity, the flux is easily computed, 

_ Txjiso 

" ^ 4^’ 

where D is the luminosity distance. 


2.3. Radio Lightcurve 

In this section, we first reproduce the lightcurve at 1.4 GHz of Sw J1644 and then we 
generalize it to events at different redshifts and with different stellar and black hole masses. 
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The radio emission is synchrotron emission and the low energy spectrum can be 
described with the following broken power-law 


F,, = FJu. 


X 


flGranot fc Sari 


Z/ «2 
1 + — 
I'm 


V 

-1/S2 


—2si 


-si/3 


-1/si 


( 8 ) 


200211 where z/a < are respectively the absorption and peak frequency, 


Si, S 2 are smoothing factors and the electron power-law index has been assumed to be 2.5. 


Berger et al.l (120121) measure the flux Fj,(z/a) = Tl/,sio(z^a,sw), and characteristic 
frequencies z^a = r'a.sw and z/m = z/m sw, in several sna pshots that cove r the e volution of the 


lightcurve up to ~ 220 days after the trigger. Later, 


Zauderer et ah 


(2013) extended the 


period of the radio monitoring up to ~ 600 days. The hrst observation is at ~ 5 days after 
the detection in the X-ray band. Therefore, the radio emission is observed after a delay 
r ~ 8 days from the intrinsic beginni ng of the event. Finally, note that the radio data 


monitoring occurs up to ~ 600 days (IBerger et ah 


2012 


Zauderer et ah 


2013 1. while the 


X-ray emission has been observed up to ~ 500 days. This mismatch, however, is not a 
problem, since we are interested in modelling the lightcurves only up to one year after the 
explosion, when is already too dim to be detected by an X-ray survey in most cases. 

Using the available data and eqJHl we can therefore model the temporal evolution of 
the flux at any radio frequency. In Figure 121 we show the lightcurve of Sw J1644 at 1.4 
GHz, and its comparison with data. A smooth temporal behavior has been obtained by 
linearly interpolating the flux between data points. 

We now need to generalize our prototypical lightcurve to a generic TDE. The main 
uncertainty is how the flux scales with black hole and stellar masses. A hrst possibility is 


to describe the jet evolution with a Blandford Mckee (thereaf 


usually adopted for y-ray burst afterglows (e.g. 


Metzger et ah 


er “BM model”) so 


2012 


Berger et al. 


ution. 


20121 ). 
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Frequencies below 5 GHz are in the self-absorbed 
whole observed duration of the event (see £g.3 in 


part of the synchr otron spectrum, for the 


Berger et al 


2012I) . The observed specihc 


luminosity in this regime [v < is given by the Raleigh Jeans part of the Black 

Body spectrum B{h'/6)6^ oc (see eq. [8]), with a kinetic temperature given by 

3khT = meC^7min, where the minimum Lorentz factor for the shocked accelerated electrons 
is 7min r. Therefore the specihc radio luminosity is 


oc B{u/6)6%r9jf oc (rr)2, 


( 9 ) 


where (r6*j)^ is the emitting area, and we are assuming T"^ ^ 9y In the blast wave modelling 
of J1644, the external medium swept up by the jet is better described by a power-law 
densi ty decay that goes as r“^, rather than a constant density environment fjZauderer et ah 
2011). This implies T oc and r oc where F) ~ Lyptmin oc m* is the total jet energy. 
Therefore eq. [9] becomes, oc (Lj^pfmin)^'^^ oc where there is no dependence on the 

black hole mass, but only on the stellar mass. 


The simple bias 


the radio spectrum flBerger et al 


wave solution, h owever, does not describe the whole evolution of 


20121) . Therefore, we also consider a simpler approach. 


In line with our treatment of the X-ray hux, we may assume that the radio luminosity is 
proportional to the jet peak luminosity oc Lj,p oc rather than to its total 

energy, (see the X-ray analogous, eqEl which bears the same mass dependencies). As an 


extra motivation, this pr escription may 


arrested” jet model (e.g. 


Naravan et al. 


3e ius tihed in the context of the “magnetically 


20031) . We will call this prescription “the Mass 


Dependent Luminosity” model (thereafter MDL model). 

The scaling of the peak hux for sources at diherent redshift, with diherent black hole 
and stellar masses (but at the same observed time r from the beginning of the event) would 
be 

( 10 ) 


Ei/{ysiiT^ Fj/,sji;(^^a,sw) "^sw) ^ ^*,1 ^ 


1 + Zs 


D 
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for the BM solution and 


P'uiya,^ T~') Fv^swiysifi'xi T~sw') ^ 1 


1 + Z 
1 + ^sw 




D 


( 11 ) 


for our second approach. The equivalent delay at which we need to calculate the flux of Sw 

iftec0 ^ 


J1644 is r„s = rl^ 


The characteristic frequencies need to be redshifteqj according to 


IZAT = l/a,.sw r, 


1 + ^ 


a,sw V' sw / I H , I ? 

1 + 


-1 


and 


J ^m,swVswJ 


1 + ^ 
1 + ^sw 


-1 


In all cases, the flux Tl/,siu(r'a,sw, hiw) and the characteristic frequencies at any time Ts„ are 
obtained by linearly interpolating the available data. For Tsw < 8 day we extrapolate the 
radio light curve to earlier epochs. 


3. Black hole mass functions 


The mass distribution of black holes as a function of redshift is an essen tial ingredien t 


to calculate TDE rates. Since black holes grow mainly by efficient accretion fjSoltan 


1982h . 


one can calculate these functions using the mass continuity equation, given a radiation 
oution of Eddington ratios. In this paper, we use the results from 


efficiencv and 

distrib 

Shankar et ah 

(2013) 


(120131 ). In particular, we consider the two accretion models which yield the 


^Formally, one would need to consider the transformation due to different Doppler factors 
between jets. However, we here assume that all jets have approximately the same Lorentz 
factor F and viewing angle of nearly ^ 1/F. The latter is because the viewing angle 
probability (oc do, between 0 < do < is the highest at F“F 
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largest and the lowest black hole comoving number density (j){M,z), and are still consistent 
with the quasar bolometric lum inosity functions an d the local black hole mass function 


(models labeled G and G{z) in 


Shankar et ah 


20131) . In this way, we can estimate the 


uncertainty due to the black hole mass distribution of our expected TDE rates. In Figure 
[3] upper panel, we show the mass distribution functions and their uncertainty strips as a 
function of redshift, for M = 10®Mq and M = 10®Mq black holes. In Figure [21 instead, we 
show the “intrinsic” TDF rate as a function of redshift, 

/•Mnjax 


R{z) = / <P{M,z)V{z)NtaJM, 

^min 


( 12 ) 


where we denote with V{z) the comoving cosmological volume. iVtde = 10 ^ yr ^ is our 


hduc ial TDE rate per galaxy: t his value is in 


201lh . 


2013) and observational claims flDonlev et ah 


le range of theoretical expectations fiMerritt 


2002 


Gezari et ah 


2001 


Van Velzen et ah 


The minimum black hole mass (here and thereafter in our calculations) is 
= lO^Mo, as just a few SMBHs have been observed with a lower massj^ 


4. Monte Carlo calculations 


Assuming the X-ray and radio modelling described in §2, we perform Monte Carlo 
simulations (MCs) to derive the number of jetted TDEs to be detected per year, for given 
flux limit and sky coverage. 


Beside the BH mass, the main ingredients of our MCs are the trigger lag time, r, and 


3The 


recen 


Maksym et al. 


disco very of TDEs in dwarf galaxies ( ^ IO^Mq) flDonato et ah! I20l4 


2014alJbl ) seems particularly promising in overcoming this limit and use TDEs 


to hnd lower mass BHs 
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the mass of the disrupted stars, m*. The former is randomly extracted from a uniform 
distribu tion b etween 0 and 1 yiQ. The latter follows a Kroupa Initial Mass Function, (IMF 
200111 ■ 


Kroupa 


/(m) IX < 


m 

0.01 < m* < 0.08, 

— 1 

m , 

0.08 < m* < 0.5, 

^-2.3, 

m* > 0.5. 


(13) 


In fact, for each black hole mass, the minimum stellar mass is set by the requirement that 
the tidal radius should be greater than the last stable orbit (we assume a non-spinning BH). 
This requirement implies that = max[0.01, O.OdSMg]. Note that for M = IO^Mq, the 

minimum mass is rn^min = 4.5M0. Therefore, events associated with high BH masses are 
suppressed in numbers by the steepness of the IMF, as only 0.4% of all stars have m* > 8. 
However, they are in average brighter, because the average m* is larger. 

In our simulation, we start by considering the intrinsic rate R{z) (eq. [T^ properly 
modihed by accounting for the relativistic beaming, which results in a reduction by a factor 
of 27rF“^/(47r) = 1/(2F^): this is the fraction of solid angle subtended by the emission, 
when considering a two sided jet. Our hducial value for the let Lorentz factor is F = 2, as 


inferred by radio observations (F 2, 


Zauderer et ah 


2013; 


Berger et ah 


2012). If the jet 


decelerates, this value has to be intended as an average one, over the observation period. 
However, we note that this is a geometrical scaling factor and our results may be easily 
re-scaled by assuming different values of the jet bulk Lorentz factor. In addition, R{z) is 
scaled for the fraction of the sky surveyed by the assumed instrument. In the calculation 
of R{z) we have adopted both G and G{z) models in order to account for the systematic 


■^we do not use longer time lags because any extrapolation beyond the currently avail¬ 
able radio data would make our estimates more mo del-dependent, since there is no hydro- 
dynamical model that can reproduce the whole radio behavior of J1644. 
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uncertainties in the mass function modellings. The number of trials in MCs is properly 
hxed by requiring a high statistics level in each mass and redshift bin (typically > 10^). 


5. Results 

In this section, we hrst validate separately our emission models for X-ray and radio 
light curves, by comparing our predicted rates with current instruments and survey results. 
In fact, we hnd that current data do not put strong constraints on our modelling, as we will 
explain in the following. Future data have instead a greater potential. In the SKA era, we 
propose that a strategy where radio will be triggering X-ray facilities can allow us not only 
to detect but also to identify and investigate jetted TDEs in a multi-wavelength fashion. In 
the following, if not otherwise mentioned, our results are derived adopting F = 2. 


5.1. Comparison with current surveys 


5.1.1. Hard X-rays 


So far, only two jetted TDE candidates have been detected by BAT, implying an 
observed rate of ~ 0.3 yr“^. 

Since BAT is not operating in survey mode, it is not straightforward to compare 
observations with our predictions, i.e. it is difficult to chose sky coverage and detection 
limit, because they are not univocally determined. The two TDE candidates were detected 
in two different modes: Sw J1644 was triggered onboard, while Sw J205 8 was discovered by 


stacking 4-day integration images flKrimm et al 


2011 


Cenko et al. 


2012 ). In both modes, it 


is hard to dehne a survey flux limit, the key ingredient of our MCs. Indeed, Swift has over 
500 onboard trigger criteria in different modes which makes the use of a flux limit survey a 
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rather simplified approach. The same applies to possible discoveries of fainter TDEs w ith 


longer integration times, by applying the image mosaics techniqne flKrimm et al 


201811 ■ 


These have to be promptly followed-np by XRT for their identification: monitor the soft 
X-ray emission and then measnre the characteristic temporal slope of TDEs. In this way, 
a farther efficiency acconnting for any reason preventing XRT to monitor the event has 
to be inclnded in onr rate calcnlations (e.g., the stochastic natnre of the Swift pointing 
plan, the target visibility and the mission schednle; Krimm, private commnnication). Snch 
an efficiency is hard to qnantify and any assnmption wonld be arbitrary and wonld bias 
onr discnssion on the comparison between the predicted and observed rates. In addition 
to that, onr soft X-ray modelling assnmes a total disrnption of the star (see §2), while 
the Sw J2058 emission seems to be consistent with a partial disrnption. For both of these 
reasons, we focns on detections triggered onboard, althongh with dne caveats. In fact, 
any relia ble prediction ba sed on on-board triggers wonld reqnire complex simnlations as 


done by flLien et al. 


2014J ) for GRB rates. We therefore set to achieve a less ambitions 
aim at predicting indicative rates, which shonld be considered most likely as nppe r limits. 


Specifically, we adopt the BAT daily sky coverage reported in flKrimm et al. 


2013 1 and fix 


a nniqne “survey” flux limit to be consistent with the detection of the Sw J1644. We detail 
the procedure in the following. 


Sw J1644 was detected with an on-board BAT image trigger fICnmmings et al. 

s 

band, which is consistent with the faint tail o 


201111 


In this trigger mode, we assume a flux limit of 2.5 x 10 ® erg cm ^ s ^ in the 15 — 


and the detec 


flKrinim et al. 


ion of Sw J1644 flBurrows et al. 


time (iLien et al. 


the observed GRB rate (iLien et al. 


50 keV 


20i4 


201111 . We adopt a daily sky coverage of 


2 OI 3 II a nd apply an efficiency of ~ 90%, for the fraction of the BAT survey 


201411 resulting from trigger deadtimes (e.g. due the passage through the 


South Atlantic Anomaly). 
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The detection rate is estimated by performing a large set o 
sections 12.21 and HI In eqEl we use a radiation efficiency of ~ 0.3 


MCs as described in 


Burrows et al 


2nill i. For 


each event, we compare the flux at the trigger time r with our flux threshold. We obtain 
a TDE rate of ~ 10 — 20 events yr“^ (see Table [T]). The rate distribution with redshift 


extends up to Zr, 


O. 341 I and peaks at z ^ 0.2. The peak value ranges between 1 — 5 


yr 


-1 


(see Tabled]). The peak of the corresponding BH mass distribution is at IO^Mq and 


contains 23% of all events. 


Given our predicted mass and 2 ; distributions for the observed TDEs, an event like Sw 
J1644 has a chance probability which is a factor of ^ 10 lower than that of an event at the 
peak rate. Therefore, it is not an unlikely event, but a lower redshift object would have 
had a higher probability. At this point, it is unclear to us if this result is more due to our 
simplihed treatment of the BAT trigger, to our assumption of a constant jet luminosity for 
a given BH and stellar mass. Both are very likely to have a role. But since we can not trust 
at this level our trigger modelling and the paucity of detected events does not constrain a 
possible luminosity function, we do not attempt here to modify our X-ray model to ht this 
observed distribution. When more events will be identihed, our procedure can be rehned to 
account for a TDE variety. A comparison with future, easier to model, surveys (see Sec. 
15.2. ip will also help constructing a more robust emission model. 

Interesting, these rates are actually up to two order of magnitude higher than that 
(0.3 yr“^) derived from BAT observations, but a key role is played by the low value of T 
considered. We will elaborate on this point in section 16.11 

®Here and in the following, we dehne Z m»^ as the redshift at which the expected rate is 
0.5 yr“^. 
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5.1.2. Radio surveys at 1.4 GHz 


We compare our predictions w ith constraints on the jetted TDE rate derived from 


current radio surveys ([Bower 


201ll) . In the following radio estimates, we will require a 5-cr 


flux limit to claim detection. 

We hrst consider the combined catalogs of VLA First and NVSS at 1.4 GHz. The 
combined sky coverage is 0.19 sr with a flux limit of 6 mJy at 1.4 GHz. The analysis of 
these catalogs didn’t yield any TDE candidate. 

To derive our predictions, we adopt the radio modelling described in section 2.3, and 
for each event (i.e. for each set of r, black hole and stellar mass, and redshift), we calculate 
the average flux over a period of one day from the trigger. This is compared with a 6 mJy 
flux threshold. Rescaling our all sky results for the catalogue sky coverage, we obtain an 


observed rate t 
consistent with 


lat even in t 


re most favorable case 


Bower! (120111 ) and 


Frail et al 


0.3 yr ^ using BM model, eq JlU]) is 


f|2012h results. To strengthen this conclusion. 


we note that our assumption of a 6 mJy threshold per day combined with a sky coverage of 
0.19 sr may be considered already rather optimistic, since both values are referred to a 1 yr 
single epoch. 


In the near future, the VLA Stripe 82 survey may constrain jetted TDE models than 


cs 


to th e improved sensitivity (50 /iJy rms) at 1.4 GHz over a FoV of 90 deg^ ([Hodee et al 


20131 ). By assuming a 5-a threshold of 0.25 mJy, our modelling predicts a number of a few 


objects to be detected per year. Signihcant advances in TDE detections are expected to 
come from on-going wide radio surveys at both low (see e.g. MWA and LOFAR) and high 
radio frequencies (e.g. ASKAP and VLASS). Since our radio modelling was constrained 
by observations at higher (> 1.4 GHz) radio frequencies (as discussed in §2.3), we focus 
here on ASKAP and VLASS. The Variable and Slow Transient (VAST) project on ASKAP 
envisages a sky coverage of 10^ deg^ reaching a sensitivity of 0.5 mJy rms (VAST wide) 
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with a daily cadence flMurphv et al. 


20131) . Our predictions for VAST are in the ran ge of 


a few up to ^ 14 TDE yr consistent with expectations from 


comparison, 


Frail et al. 


Murphy et ah 


(120131). For 


fl2012[) obtain a value of ~ 82 yr 


-1 


by considering long er integration 


2013 


, VLASS), we 


time ten days). In the case of the VLA Sky Survey (iHaHinan et al. 
consider a sky coverage of 10^ deg^ with ~ 3 week cadence at a sensitivity of 0.1 mJy rms. 
This set up should give a number between 2 and 6 objects to be detected per year. All sky 
VLASS is also foreseen and will clearly provide a larger number of TDEs, but we focus on 
the previous strategy because the multi-epoch survey could provide alerts for follow-up at 
higher energies, with a prompt identihcation of the transient. 


5.2. Future instruments 

Currently, the only two jetted TDE candidates were discovered in X-rays, where the 
characteristic decay slope has been observed. Therefore, we hrst discuss the discovery 
potential of future X-ray surveys. We then predict the expected rate of TDEs for the SKA 
1.4 GHz wide survey. Finally, we derive the properties and rate of TDEs that can be 
detected in radio with SKA and subsequently identihed in X-rays. 


5.2.1. Future X-ray surveys 

The rate estimates provided in this section are based on a unique observing strategy 
aimed at detecting and providing a hrst identihcation of the transient as a TDE. We 
assume that a given fraction of the sky is covered in 1 day at a hux threshold dehned by 
the requirement to follow the typical TDE decay over 4 lightcurve bins, each with S/N > 5 
This is obtained by starting from the 5-0" hux limit of each survey, then tracing back the 
f-5/3 order to obtain the hux over the 4 time bins and then compute the average 
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flux over that period. This average flux defines the identification flux threshold. We will 
give values for both T = 2 and T = 20 and we will justify this choice and elaborate on the 
comparison in section O 


The all sky survey mission eRosita (iMerloni et al.l 120121) is expected to detect jetted 
TDEs, in its “hard” X-ray band (2 — 10 keV). We apply our methodology to the eRosita 
survey, properly re-scaling the sky coverage achieved in a 6-month scan to 1 day. We 
derive the identification flux threshold for our observing strategy from the 2 — 10 keV 5-cr 


sensitivity of ~ 10 erg cm ^ s corresponding to ~ 250 s exposure flMerloni et ah 


201211 


as foreseen for each point in the sky. We calculate the corresponding un-absorbed flux 
and then we extrapolate it in the energy range 1 — 10 keV (used in our X-ray modelling). 
We predict a maximum of ~ 15 TDE per year to be detected up to .2 ~ 2.5, although 
^max ~ 0.4. The peak rate is between 0.15 and 0.5 yr“^ at 2 ; ~ 0.4 and beyond 2 ; ~ 2, the 
rate is < 5 x 10“^ yr~^- If a larger value of T = 20 is considered, the rate decrease by two 
orders of magnitude (see Tabled]) with a maximum total rate of ~ 0.15 yr“^ and peak 
rate of only 4 x 10“^ We therefore predict both higher (T = 2) and lower (T = 20) 


rates than those previously published by 


KhabibuHin et ah 


(12014 1 object to be detected 


per 6-month long scan), but we dehnitively reach a much lower redshift (zmax = 0.4 vs 
their 2 ; = 4.5). The same authors provide an upper limit of ~ 150 events per scan by 
considering the number of jetted TDEs to be a l/5th of their “soft” TDE sample 1000). 


This fraction is based on results from the Rosat X-ray survey flDonlev et ah 
clearly consistent with their estimate. 


2 OO 2 I) . We are 


The Wide Field Monitor (WFM) aboard of a LOFT-like ( Feroci et al.ll2012l ) mission 


will also have the capability to trigger jetted TDEs by surveying l/3rd of the sky with a 
5-a 1 day sensitivity of ~ 8 — 9 x 10“^^ erg cm“^ s“^ (a few mCrab) in 2 — 50 keV energy 
band. We estimate tens of objects per year up to ~ 0.6. The peak rate is ~ 6 yr“^ at 
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2 ; = 0.2. These numbers imply a total rate of 0.7 yr ^ for T = 20 with a peak rate of only 

6 X 10“^ yr“^- 


Finally, we consider Einstein Probe, which is expected to monitor 1/2 of the entire 
sky in the energy range 0.5 — 4 keV with a 5-cr sensitivity of ~ 10“^^erg cm“^ s“^ in each 
point (1 ks exposed) of the sky (EP, W. Yuan private communication). In this case, MCs 
were adapted in order to extend our X-ray modelling to this energy range. This requires 


to hrst estimate the un-absorbed f 


intrinsic absorption (jBnrrows et al 


ux lim it by accounting for both the Galactic and the 


20 111 ) and then calculate a proper radiation efficiency 


by extrapolating from the value inferred in 1-10 keV. We estimate a number between 
~ 90 — 240 yr“^ to be detected below ~ 1 with a peak rate of ~ 15 yr“^ at z = 0.3. In 
the case of P = 20, a few objects are expected to be detected per year, with a peak rate of 
~ 0.2 yr“^. A summary of the actual numbers can be found in Table [H 

Inspecting the trigger time distributions (see top panel in Figure [5]), we hnd that up to 
Zmax ~ 1 objects are detected with almost equal probability at any delay from the explosion. 


These X-ray survey rates have been obtained under the assumption of a reasonable 
observing strategy. A larger sample extending up to higher redshift can be obtained if longer 
integration times are considered, but these predictions are affected by several parameters 
like the trade-off between sky coverage and sensitivity. In this respect, our approach has to 
be considered conservative. 


5.2.2. SKA as TDEs hunter 


Presently, the most ambit ious and revolutionary p roject in radio astronomy is the 
Square Kilometer Array (SKA 
in survey mode (SKAl-Survey, 


Carilli & Rawlines 

2004 

) planned to operate in 2020. SKA, 

Dewdnev et al. 

(2 

013 

)) 

is able achieve a half sky coverage 
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(20,000 deg^) with 

a 2-d< 

2014; 

Feretti et ah 

— ^ 
r—1 

0 

CM 


ay cadence at a 5 — a flux limit of 90 /iJy flDonnarnmma et a1 


2014J) . These unprecedented sky coverage and sensitivity make SKA an 


optimal radio transient hunter. 


Differently from X-ray searches, in radio, we cannot have a first identification based 
on the lightcurve, since the 1.4 GHz radio emission of a TDE is not particularly different 
from those of other radio transients (e.g. GRB, blazars). Therefore, we consider a different 
strategy. In our MC simulations, we directly assume the SKA 5 — a flux limit in order 
to claim the detection of a transient event. The identification strategy will fully rely on 
the multi-frequency follow-up of the trigger event as it will be discussed at the end of this 
section. 


We calculate the predicted average flux over 2 days from the trigger and then we 
compare it with the SKA flux limit. The results are shown in Figure HJ The upper panels 
are derived using the BM model feq JTOl) for the radio lightcurve modelling, while the lower 
panels use the MDL model (eq ITTll . There, we show the distribution of the TDE rate as a 
function of 2 ; (right panels) and their BH mass distribution (left panels) for the two BH mass 
functions described in §4 (black lines). The yellow lines show the subclass of events with 
BH masses lower than 10^ Mq. Both radio models produce redshift distributions peaking 
around 2 ; ps 0.4, regardless of the BH mass function. The peak rates are roughly between 
6 and 40 yr“^ (see also Table [T]). Events with BH mass lower than 10^ Mq dominate the 
distributions at all redshifts in the MDL model, while this only happens at z < 0.4 in the 
BM model. One marked difference between the two radio lightcurve modellings is the BH 
mass distribution of the detected events: while BH with masses between 10^ — IO^Mq are 
equally probable in the BM model (because the flux is BH mass independent, eq JTOj) . BHs 
with mass < 3 x 10® completely dominate the observed sample in the MDL model. As a 
consequence, BM model distribution extends to higher redshifts (z m^^ ~ 2.5 vs z^ax ~ 1-7), 
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because BH masses larger than 10’^ Mq interact with higher mass stars (higher and 

produce intrinsically brighter flares. This will allow us to study TDEs close to the peak of 
cosmic star formation. 


In Table [H we also report the total rates obtained by integrating these distributions 
in .2 and Mbh- We obtain yearly rates of the order of a few to several h undreds. These 


resul ts are not consistent with those that can be derived by using eq.4 in fjVan Velzen et ah 


20131) : inserting our SKA survey parameters, we obtain thousands of events per year. This 


discrepancy is due to our inclusion of the stellar mass dependence, that modulates the TDE 
luminosity for a given BH mass: the lower m*, the dimmer the event. In the assumption of 
a Kroupa IMF, the bulk of the events are caused by the disruption of stars with m* < 1, 
increasing the number of flares that are too dim to be detected. 


With hundreds of events per year, SKA could be able to detect more TDEs than any 
currently planned X-ray survey. On the other hand, while X-ray surveys can catch the 
events soon after explosion (see EP performance in Figure [5] upper panel, for an example), 
SKA would not be able to cover the first week activity at any redshift and only at {z < 0.8) 
SKA will probe the first month (Figure O bottom panel). This result is independent on 
the assumed radio modelling. The explanation is simple: the observed radio flux at 1.4 
GHz is initially increasing, contrary to that in the X-ray band. In this regard, detections 
in these two bands are complementary. However, a word of caution here is due. As 
mentioned before, below 10 days, we have virtually any detection in radio at any 2 ;. This 
early period coincides with the rise of the radio light curve. Although we expect this gap in 
detection, the exact epoch at which it occurs depends on the detailed behavior of the light 
curve during this undetected rise. Our extrapolation at earlier times is quite steep and we 
consider 10 days as an upper limit for the initial gap in detection. 


So far we focused on detection of TDEs with SKA, that, depending on the observing 
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strategy, will only be a fraction of a noticeable sample of slow radio transients. As 
mentioned earlier, we cannot use radio properties or variability alone to distinguish a TDE 
candidate from neither a slowly variable AGN or a GRB. A possibility for identification 
that we explore below is through quick follow-ups at higher energies, particularly in X-rays. 
A Erst pre-screening of the radio candidates could be done by cross-correlating the radio 
transient positions with deep AGN catalogues, expected to be provided in the near future 
by optical surveys (e.g. LSST) or the SKA precursors (e.g. ASKAP). However, we expect 
a larger degree of contamination of the TDE sample to come from transient sources such 
as GRBs. Since, unlike GRBs, most of TDEs should have a nuclear origin, it is mandatory 
to quickly identify the host galaxy. An accurate localization of the radio transient in the 
core of galactic nuclei, helping to assess the nuclear origin, will therefore play a major role 
in the screening of the radio transient sample. This means that hrst the host galaxy has to 
be found by a rapid optical follow-up and after the brighter transients could be localized by 
SKA with a precision of ~ 100 milliarcseconqj (mas) essential to separate nuclear transients 
from other phenomena (e.g., GRB). 


For details see 


Donnarnmma et ah 


(120i4 . 


5.2.3. Combining Radio and X-rays in the SKA era 


X-ray follow-up will have a major role in the identihcation of the TDE candidate 
detected by SKA because of the possibility to detect the characteristic t~^^^ decay. A 
possible X-ray follow-up strategy aimed at identifying and then characterizing the event 
consists in a fast repointing of the transient detected by SKA. We consider a 1-day delay in 


®this can be achieved thanks to the resolution o ' 

arcsec or better of SKAl-MID ( 

Dewdnev et ah 

201s5) 


about 2 arcsec of SKAl-SUR and 0.6 
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the X-ray repointing and reqnire a set of X-ray observations spread over a few days in order 
to follow the characteristic temporal decay of the TDE. We foresee an observing strategy 
which is similar to the one adopted in the case of fntnre X-ray surveys (see section [5.2.ip : 
four observations spread over ^ days, with S/N ratio > 10 in each. A high S/N is required 
in order to characterize both the temporal and spectral behavior of the source. 

For each event in the MCs, we calculate the average X-ray flux over the 4 days after 
the repointing and compare it with the identification flux threshold derived as explained 
in section 15.2.11 with the only difference of a S/N > 10 requested in each observation. 
Practically, r in eqElhas to be the radio trigger time-lag, plus an extra delay of one day for 
repointing, and 0 < At < 3. In this way, we derive the properties of samples of TDEs which 
are hrst detected in radio and promptly followed-up in X-rays. 

In Fig. ini we show the fraction of SKA candidates that can be identihed as a 
function of the X-ray (1-10 keV) unabsorbed flux limit. A rapid X-ray follow-up will be 
able to detect a complete radio-selected sample provided that the instrument sensitivity 
is close to Eiim < 10“^^ erg cm“^ s“^ in the energy. In fact, a moderate sensitivity 
~ 10“^^ — 10“^° erg cm“^ s“^ is already enough to detect equal or a larger number of events 
than with X-ray wide sky instruments alone. It is therefore clear that a radio trigger is a 
more efficient way to build up a large X-ray sample of TDEs. Rates reported in that Figure 
assume a fast (1 day) X-ray repointing and Elim reached with an integration of ~ 4 days. 
Rates could be substantially different if longer integrations are needed to reach the same 
Fiim or in the case of longer repointing time. This is a natural consequence of the decreasing 
trend of the X-ray light curve. 

When considering an actual follow-up strategy, the values reported in Figure [6] should 
be scaled by the fraction of sky accessible to the X-ray instrument considered. In general, 
the X-ray follow-up will provide us with a sub-sample of radio triggered TDEs, dehned by 
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the target accessibility, the repointing chance of the X-ray satellite and the sensitivity of 
the instruments. Since TDEs also emit in hard X-rays, a trade-off between sensitivity, sky 
coverage and a broad energy range is foreseen. In particular, the broader is the energy range 
the better the characterization of the non-thermal process and of the jet energy budget. 


Future X-ray 


flFeroci et al. 


experiments like Athena flNandra et al. 


2013 1 and a LOFT-like mission 


2 OI 2 I) could offer a unique chance to follow-up and characterize SKA triggered 


TDFs. Moreover, if Swift were still operating in the 2020s, XRT will have a great potential 
in following-up the radio candidates. 


Athena sensitivity lies in the saturation branch of Figure [6l which implies that the 
observed rate of X-ray jetted TDEs will be crucially linked to its follow-up efficiency. 
This is mainly influenced by the Athena sky accessibility which is of the order of ~ 50% 
(Athena mission proposal), resulting in a rate of TDEs of a few hundreds, with detections 
up to Zmax ~ 2 (scc Table [T]) . 


The LAD (Large Area Detector) on board of LOFT (2-50 keV) is a collimated 
instrument with 1 degree held of view, and a background limited sensitivity of 
> 10“^^ erg cm“^ s“^ in the 2-10 keV band, for a 100 ks exposure. The LOFT pointing 
visibility will assure a sky accessibility for these targets of ~ 75%, (LOFT Yellow Book). 
The requirements of our strategy dehne a 10 erg cm ^ s ^ in the 2-10 keV band, 

which was then translated in the corresponding un-absorbed value in the 1-10 keV band 
(the energy range adopted in our modelling). Again, we assume a 1-day repointing delay. 
Figure [3, shows the expected rate of jetted TDEs for a LOFT-like mission as a function 
of redshift (right panel) and their mass distribution (left panel). The rate distributions 
are calculated under the BM model (top panels) and MDL model (bottom panels) for the 
radio modelling. In both cases, we found that the redshift distribution extends above z = 1 
(^max ~ 1-2 — 1.7) (see Tabled]), with most of the TDEs expected around ^ ^ 0.4 (right 
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panels). The peak rates are roughly between 4 and 20 events per year. In the MDL 
model, because of the mass dependence of both the radio and X-ray luminosities, ~ 25% of 
all events have BHs with masses ~ 10® Mq, and events with BH masses <10^ Mq dominate 
the redshift distribution at all epochs. In the BM model, instead, the TDE rate peaks at 
~ 10^ Mq (see left panel in Figure [7]), with lighter BHs dominating a.t z < Zpeak (yellow 
lines in Figure [71 right panel). The behavior at higher fully reflects the one observed in 
the BM radio rates (see £g. H]). In total, a LOFT-like mission should be able to detect a 
sub-sample of radio TDEs between 130 and ~ 350 yr“^. Instead, very few objects per 
year are predicted if T = 20. For these events, the mission broad energy band (2 — 50 keV) 
should enable us to put tighter constraints on the energy budget of the X-ray component, 
than possible with Athena instrument. 

The price to pay for detecting more X-ray TDFs with a follow-up strategy is illustrated 
in Figure [5] upper panel, where we compare the trigger distribution for the FP (black lines) 
and the LOFT radio triggered (blue lines) samples. Most of LOFT events are observed 
after 10 days from the beginning of the emissioijll. In particular, high redshift ~ levents 
are all a couple of months old. Direct discovery of TDFs in X-rays is thus important for 
catching the event in its very early dynamical stages, when the jet has just formed and the 
disc may still be in the (largely unconstrained) super-Eddington regime. 


6. Discussion 

The Swift/BAT discovery of Sw J1644 opened a window on a new class of X-ray and 
radio transients, which are optimal targets for future radio and X-ray surveys/instruments. 
The study of these objects allows us to investigate the formation of transient jets in 


^See discussion in Sec. 15.2.21 
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yr-i 

yr-1 

'^peak 

tdI 

■^peak 

yr-i 

d2 

■^peak 

yr-i 

^max 

Radio selected sample 







SKA BM 

226 

468 

0.3 

6 

17 

2.5 

SKA MDL 

327 

770 

0.4 

14 

40 

1.7 

LOFT-like BM 

128 (2.5) 

305 (6.5) 

0.3 

4.5 (0.05) 

13 (0.1) 

1.7 

LOFT-like MDL 

135 (1.3) 

352 (3.5) 

0.4 

8 (0.08) 

22 (0.2) 

1.2 

Athena BM 

113 (1) 

234 (2.3) 

0.3 

3 (0.03) 

8.5 (0.09) 

2 

Athena MDL 

163 (1.6) 

385 (4) 

0.4 

7 (0.07) 

20 (0.2) 

1.4 

X-ray surveys 







BAT3 

9.5 (0.095) 

26.5 (0.26) 

0.1 -0.2 

1.7 (0.02) 

4.6 (0.05) 

0.32 

eRosita 

8 (0.08) 

15 (0.15) 

0.4 

0.15 (0.001) 

0.5 (0.005) 

0.4 

Einstein Probe 

89 (0.9) 

242 (2.4) 

0.3 

5.5 (0.05) 

15 (0.2) 

1 

LOFT-like WFM 

24.5 (0.2) 

67 (0.7) 

0.2 

2.3 (0.02) 

6 (0.06) 

0.6 


Table 1: Future Radio and X-ray surveys predictions: 1st and 2nd columns are total yearly 
rate (the subscripts ^ and ^ are for the Gz and G MFs, respectively), 3rd column redshift at 
the peak rate, 4th and 5th columns maximum peak rate and 6th column maximum redshift, 
defined as the z where the rate is 0.5. BAT^: calculation for an on-board image trigger. 
X-ray and Radio expected rates are derived for F = 2. X-ray rates are also reported for 
F = 20 in parenthesis. 
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extra-galactic sources. Moreover, there is the potential to discover quiescent SMBHs in 
distant galaxies and constrain the SMBH mass function. In this section, we qualitatively 
discuss our results and what we may learn from them. Any quantitative parameter 
investigation (for instance with a Fisher Matrix technique) is beyond the scope of this 
present paper, and will be presented in a follow-up work. 


6.1. Jet efficiencies and bulk Lorentz factor 

So far, only two jetted TDEs have been detected, while the thermal candidates, related 
to the presence of an accretion disk, have been more numerous. The question then arises 
whether this is due to observational biases, highly collimated jets or to an intrinsic low 
efficiency of transient accretion disks to produce (luminous) jets. 

To try and address this question, we could compare our predictions to the Swift/BAT 
observed rate (~ 0.3 yr“^): our lower limit (~ 9 yr“^) is a factor of 30 higher. It is tempting 
— and indeed it has been done in the literature — to reconcile this discrepancy by invoking 
a jet production efficiency of a few percent, since our calculations assume that each TDE is 
accompanied by a jetO 

However, there are several reasons why this inference should not be drawn. First, as 
discussed in Section [5.l.H it is absolutely non-trivial to describe the characteristics (e.g. flux 
limit and sky coverage) of an effective Swift/BAT survey. We believe that our assumptions 
for the trigger, together with a 100% identihcation efficiency gives rates that are indicative 

®In our simplihed description here, there are only two kinds of possible events: Sw J1644 
with its own jet luminosity (i.e. a given jet energy efficiency ej) and events with no jet (i.e. 
€j very small). In reality, there must be an intrinsic distribution of ej, with a tail of low 
energy events that cannot be detected or failed to be launched at relativistic speeds. 
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of an upper limit. Second, BAT rate predictions, unlike those of other X-ray instruments 
consider here, strongly depend on the modelling of the early stage variability of the X-ray 
lightcurve (see §2). The onboard threshold we use is very close to the flux of the upper 
envelope of the lightcurve. We are therefore implicitly assuming that we can always trigger 
an event, by catching it at its maximum. However, since the flux varies by two orders of 
magnitude, our choice implies again an upper limit estimate of BAT rates. Finally, even if 
we trust our modelling of the BAT trigger and initial X-ray variability, uncertainties in the 
value of T can account for the discrepancy. So far, we have considered a bulk Lorentz factor 
of 2, since the radio measurements strongly support such a low (T ~ 2 — 5) value. However , 


hard X-ray observations are consistent with larger Lorentz factors (T < 20 


Burrows et ah 


20111 ). which will bring down our rates to the observed value (see Tab.l). The consequence 
would be that the simultaneous hard X-ray and the radio emissi ons need to come from 


different regions — as already claimed (e.g. 


Zauderer et al. 


20111) . The picture may be 


that while the radio emission is produced from further out, after the 


decelerated. X-rays probes regions much closer to the central engine (iBloom et al. 


jet has subs 


antia lly 


201ll). If 


that was true. X-ray detections and follow-ups would be further suppressed with respect to 
the expected SKA performance. 

Unlike the previous comparison with BAT results, our predictions of the radio rates 


are c onsistent with the upper limits derived using with the NVSS -|- FIRST catalog flBower 


201 ll) . for any F > 2. As a consequence, this comparison cannot provide us with further 
constraints on either F or the jet efficiency. In the next future, surveys such as VLA Stripe 


82, ASKAP and VLASS will give tigh 


sensitivity (50 /iJy rms. 


Hodge et al. 


er co nstraints on jetted TDEs thanks to the improved 


20131) of the former and the wide held of view of the 
latter two surveys. In this case, our radio modelling predicts a number of a few objects yr“^ 
(a few tens yr“^) to be detected by assuming F = 2. Comparing predictions with (positive) 
observations will thus constrain possible combinations of F and jet production efficiency. 
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As already discussed, the optical transient surveys Pan-STARRS and LSST are 
expected to make significant advances in the study of TDEs. LSST will be a real 
breakthrough in this respect, surveying 2 x lO'^ square degrees of the southern sky. 
Thousands of objects are expected to be discovered at z < 1 by catching their thermal 
light from the accretion disc or from the non-relativistic wi nd in the Super-Eddington 


phase, surveying 


Van Velzen et al 


he sam e fraction of the sky every 3 days (jStrubbe fc Quataert 


2009 


2011a|) . However, optical extinction in galactic nuclei still introduce an 


observational bias in the TDEs discovery alt 


rough less significant with respect to that 


occurring in the UV band. As suggested by flStrubbe fc Quataert 


20091) . infrared surveys 


will provide a complementary approach being the lower frequency energy range less affected 
by any source of obscuration. 


Contrary to radio and X-ray emissions, the optical and infrared light are not expected 
to be relativistically beamed nor to be connected with jet emission. These features imply 
that a comparison between optical. X-ray and radio selected samples can help constraining 
both the TDE efficiency to produce jets and the relativistic Lorentz factor. This latter, 
when an X-ray sample is available, will help assessing the jet energy efficiency ej. 


6.2. Supermassive BH masses 

To understand supermassive BH cosmic growth and their connection with the host 
galaxy, it is necessary to have a good understanding of which mass can be found in which 
galaxy and, more broadly, of the SMBH mass function as a function of redshift. 

The detection and light modelling of a TDE event is a unique way to constrain the mass 
of an otherwise quiescent BH, that is too distant to be detected by stellar dynamics. An 
attractive feature is that TDEs may occur in any type of galaxy, allowing for the detection 
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of a broader range of SMBH hosts. For the lightcurve modelling, a multi-wavelength 
approach can yield tighter constraints on the mass, since other parameters such as the jet 
energy, Lorentz factor and the stellar mass need to be simultaneously determined. 


A perhaps more direct measurement of the BH mass can come from 


very fast X-ray 


variability, as the quasi periodic oscillation (QPO) observed in Sw J1644 (iReis et al.l 120121) . 


The prospect for detection of QPOs in such events is quite favorable for both Athena and a 
LOFT-like mission. If QPOs in TDEs were associated w ith the Keplerian frequency at the 
innermost stable orbit fas discussed in 


Reis et ah 


20121) . the highest rest frame frequency 


should be of the order of 200s for a BH mass of 10® Mf:^. This QPO frequency is easily withi n 


reach of both LOFT-like and Athena instruments flFeretti et ah 


2014J: 


Nandra et ah 


201311 . 


Longer oscillations are expected for more massive BHs (oc -^/IO^/Mbh), whose detectability 
could be more complicated due to satellite orbit constraints (e.g. Earth occultation. South 
Atlantic Anomaly). However, providing that the QPO is persistent over a long period and 
the source is bright enough to remain above threshold for several cycles, a direct measure 
of such a QPO is also possible. 


An other method to constrain the mass function may be to compare our rate 
distributions with future SKA triggered observations. As shown by Figure |3] upper 
panel, there are still uncertainties in the BH mass function, which in turn affect our rate 
predictions (see Figure |3] lower panel. Figure [Hand Figure [7|). 


7. Conclusions 

We have investigated the best strategies to increase the sample of the new class of 
TDEs, which was recently discovered by BAT. These events emitted non-thermal emission 
in X-ray and radio bands, probing a relativistic jet. Given the lack of statistics and 
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of a solid theoretical framework for their non-thermal emission, we adopted a rather 
phenomenological approach to model their lightcurve. We £t the behavior of the best 
stndied candidate, Sw J1644, in both radio (1.4 GHz) and X-rays (1-10 keV), and we 
nsed the classical theory of TDEs to rescale the emission for different black hole and star 
masses. In the radio band, we also considered, in alternative, the blast wave model, usnally 
adopted for GRBs. We then nsed a Monte Garlo code to compnte their expected rate as 
a fnnction of redshift and black hole mass. We considered both current and future radio 
and X-ray surveys/instruments. Since the characteristic temporal decay of a TDE event 
can be observed in X-ray, an identihcation is claimed only when the X-ray emission can 
be sampled in at least 4 lightcurve bins with high signal to noise ratio, S/N > 5. When 
the TDE is detected in radio, we investigated a follow-up strategy for identification which 
required X-ray detectors to sample the lightcurve with the almost the same requirements as 
above (but with a S/N > 10 ). To concretely explore future possibilities, we investigated 
in particular the expected performance of eRosita, Einstein Probe, Athena, a LOFT-like 
mission and SKA operating in survey mode (SKAl-SUR). 

Our major hndings can be summarized as follows: 


results from current instruments (such as BAT and NVSS -|- FIRST catalogues) do 
not provide constraints on jet parameters or the jet production efficiency; 


However, to reconcile BAT predictions wit h observations a T 


consistently with hard X-ray observations flBurrows et ah 


20 may be adopted. 


20111) . If this were true. 


X-ray and radio e missions should come 


on different bases fjZauderer et ah 


rom two different regions, as already suggested 


20111). The predicted X-ray rates would also be 


suppressed by (2/20)^ with respect to those in the radio band. 




In the near future, VLA Stripe 82 survey, VLASS and ASKAP-VAST may provide 
from a few to ten events yr“^, putting some constraints on possible combinations of 
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bulk Lorentz T and jet production efficiency; 


Hundreds (F = 2) of Sw J1644-like objects per yr are expected to be within reach of 
SKAl-SUR at 1.4 GHz. They can probe the distant Universe up to z ~ 2.5. These 
results differ froin prev ious, more optimistic, predictions of thousands yr“^ (for T = 2 


Van Velzen et af 


20111 ) 


• Future X-ray surveys will provide a more modest sample, between several (eRosita) 
to a maximum of ~ 240 (EP) jetted events per year. With a highly collimated jet, 
with T = 20, these numbers drop to a maximum of a few. 


• X-ray detections can be substantially enhanced, if a prompt follow-up of SKA 
candidate is adopted with an instrument with flux limit A. 10“^^ erg cm“^ s“^ in the 
1-10 keV band over 4-day timescale. With that flux limit each SKA triggered event 
can have in principle an X-ray counterpart (see Fig|6]). A suppression factor should 
be adopted if the X-ray emitting region would be moving with a larger Lorentz factor. 

• The sample of SKA preselected X-ray events can extend up to redshift ~ 2 for a X-ray 
instrument such as Athena and the LAD on board of a LOFT-like mission. Instead, 
eRosita, the WFM on LOFT and EP samples will probe a redshift range only up to 
z 


• Despite the several advantages of a radio trigger, direct X-ray detections are the only 
way to study the early stages (< 10 day) of the flare (see Figure |5]). 


Once TDE samples in different bands have been built up, the synergy between radio. 
X-rays and optical can in principle constrain important physical quantities such as the 
jet luminosity, bulk Lorentz factor, the jet production efficiency and the black hole mass 
function. These Endings will inform theories of jet and disc formation from sudden accretion 
events and, on the other hand, of SMBH cosmological evolution. 
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Fig. 1.— The X-ray (0.3-10 keV) lightcurve of J1644 as a function of time from the X-ray 
trigger: data (absorbed flux, circle marks taken from Publicy available XRT lightcurves) 
versus our modelling (solid line). After a few days, the temporal decay approaches 
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Fig. 2.— The radio (1.4 GHz) light curve of SwJ 1644 as a function of ti me from the 


radio trigger (5 days after 


2012 


Zauderer et af 


he X-ray hrst detection): data (circle marks) from (iBerger et ah 


20131 ) versus our modelling (solid lin e). While our mo delling well 


reproduce higher radio frequencies lightcurves (see £g. 1 in flBerger et al. 


2012)), it slightly 


underpredicts the 1.4 GHz one. In this respect our flux modelling is conservative. 
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Gz vs G models 




Fig. 3.— Upper panel; BH number density as a function of .d for 10® (black shaded area) and 
10® Mq (blue shaded area). Lower panel; intrinsic rate of TDEs as a function of redshift. A 
rate of 10“® yr“^ per galaxy is assumed. Most of the events are expected below z ^ 2. 
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Fig. 4.— Rate of events predicted for SKA in wide survey mode at 1.4 GHz as a function 
of redshift (right panels) and their distribution as a function of BH mass (left panels), for 
two different black hole distribution functions (black solid line: G model black dotted line: 
Gz) model. Rates associated to events with BH masses lower than lO'^ Mq are also shown 
(yellow lines). Upper panels: BM model for the jet evolution; lower panels, MDL model. 
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Fig. 5.— Cumulative distributions of delays in detecting the TDE from the explosion time, 
for different redshifts. Top panel: EP (black lines) and LAD follow-ups of radio triggered 
TDEs (blue lines). The different line styles are for 2 ; = 0.1, 0.2, 0.37, 0.8 from right to left for 
EP and viceversa for LAD. Bottom panel: the same as above but for SKA BM model and 
z = 0.1, 0.2, 0.37, 0.8,1.5, 3 from left to right). 
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Fig. 6.— Fraction of X-ray TDEs that can be identified, following up a SKA trigger, as a 
function of flux limit (unabsorbed flux). The blue and black shaded area are obtained with 
eqfTOl and eqUU radio modellings, respectively. The shaded areas reflect uncertainties in BH 
mass functions. The figure shows that a X-ray instrument with a flux limit of ^ 10“^^ in 


cgs units, can in principle identify any radio detected TDE. See text for details. 
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Fig. 7.— As Figure m but for radio triggered events repointed by a LOFT-like mission. 































































